Magnetic properties of nanoparticles in the Bethe-Peierls approximation 
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In this work we present a new method to calculate the classical magnetic properties of single- 
domain nanoparticles. Based on the Bethe-Peierls (pair) approximation, we developed a simple 
system of equations for the classical magnetization of spins at any position within the nanoparti- 
cle. Nearest neighbor pair correlations are treated exactly for Ising spins, and the method can be 
generalized for various lattice symmetries. The master equation is solved for the Glauber dynamics 
(single-spin-flip) in order to obtain the time evolution of the magnetization. The capabilities of the 
model are demonstrated through nontrivial calculations of hysteresis loops as well as field cooling 
(FC) and zero field cooling (ZFC) magnetization curves of heterogeneous non-interacting nanopar- 
ticles. The present method automatically incorporates the temperature and could be adapted to 
describe an ensemble of interacting nanoparticles. 
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I. INTRODUCTION 

The properties of magnetic nanoparticles and fine 
particle assemblies have been attracting the interest of 
physicists for many decades. Below a critical size such 
nanoscopic magnetic systems develop new properties that 
cannot be found in macroscopic systems, where the ma- 
terial is divided into domains in order to decrease the 
magnetostatic energy. Stoner and Wohlfarth |Q laid the 
theoretical basis, which are widely used nowadays, to de- 
scribe such low dimensional systems and predicted new 
effects like the superparamagnetism. According to the 
Stoner- Wohlfarth's model, the strong exchange interac- 
tion aligns all the particle spins so that the magnetization 
inside the nanoparticle is assumed to be a uniform field, 
which rotates coherently due to the presence of an exter- 
nal magnetic field. Despite being a suitable assumption 
for homogeneous systems, it is not appropriate to de- 
scribe heterogeneous particles, for instance, formed by 
regions of localized moments having distinct magnetic 
character (ferromagnetic, ferrimagnetic, and antiferro- 
magnetic). For such heterogeneous systems the exchange 
energy plays a central role in determining the magneti- 
zation of the nanoparticle (2|-Q . Unfortunately, the non- 
trivial interaction among the spins is responsible for a 
highly correlated behavior that is very difficult to de- 
scribe in terms of analytical expressions or exact numer- 
ical solutions, for a general case. In addition, the bound- 
ary conditions imposed by the finite size of the nanopar- 
ticles and surface disorder complicate otherwise simpler 
solutions for infinite systems. 

In most experiments the magnetization due to an en- 
semble of nanoparticles is the measured parameter. In 
this case the long range dipolar interaction is dominant; 
it produces a collective energy barrier for the magnetiza- 
tion reversal of a nanoparticle and, consequently, sets a 
new time scale for irreversible processes Q . The magneti- 
zation of independent nanoparticles being experimentally 
accessible in the limit of very dilute particle ensembles. 
However, a series of recent experiments, using a micro- 



SQUID technique, have demonstrated the possibility of 
measuring both the magnetization reversal mechanism 
and dynamics of individual homogeneous nanoparticles 
(10-20 nm in diameter) and cobalt nanoclusters (3- 
5 nm in diameter) j^]. The formalism presented in this 
paper should provide a theoretical foundation for exper- 
iments made on single (or independent) heterogeneous 
nanoparticles. 

Magnetic systems are well known for developing long 
range interactions, specially in the continuous phase tran- 
sitions, where the correlation length is infinite. Because 
of that the calculation of critical exponents require so- 
phisticated theoretical formalisms such as renormaliza- 
tion group techniques or high temperature series expan- 
sions, or numerical methods such as the Monte Carlo. 
On the other hand, away from the critical point, where 
the correlation length decays exponentially, a mean field 
theory may provide a good description of the behavior of 
the system. Nevertheless, the plain mean field approach 
corresponds to one in which there is no correlation be- 
tween any pair of spins and a given spin interacts with 
the average magnetization of its neighbors or, equiva- 
lently, with an effective mean field produced by all the 
spins in the system. This simplified approach might fail 
however for nanoparticles having a core with magnetic 
properties that differ from those of the surface, as in the 
instance of a ferromagnetic particle coated by a ferrimag- 
netic or antiferromagnetic layer caused by oxidation 
processes, or the case of a disordered surface as well [||. 
For these systems the correlations between neighbors are 
important to establish boundaries between regions of dif- 
ferent magnetic properties. 

In this work we propose a new method to describe the 
classical magnetic properties of non-interacting nanopar- 
ticles. The method takes a further step from the stan- 
dard mean field approach and treats the nearest neighbor 
pair correlations between spins within the particle explic- 
itly, therefore being capable of describing heterogeneous 
systems. In spite of being extensively used to calculate 
the magnetic properties of infinite systems, the Bethe- 
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Peierls (also called pair) approximation has not yet been 
used to calculate the magnetic properties of nanoparti- 
cles. A somewhat similar approach was used to describe 
the magnetization of Fe quasi- ID clusters ||^, in view of 
the fact that the results obtained by the Bethe-Peierls 
approximation agree with the Ising solution for ID sys- 
tems ^] . In the following the details of the model and 
the formalism are presented. Then, in the remainder of 
the paper, it is shown that the model accounts for many 
of the properties of magnetic nanoparticle systems. 



II. THEORY 

At first, consider a two-dimensional (2D) array of Ising 
spins in an hexagonal lattice, as depicted in Figure 1 for 
a structure made up of 3 shells. A similar model sys- 
tem could be used to describe the properties of single- 
crystalline disks with hexagonal contour, as the 10-20 nm 
BaFeCoTiO nanoparticles whose magnetization reversal 
has been measured ||]. We note that the formalism ap- 
plies equally well for any 2D lattice and can be extended 
to tridimensional (3D) lattices as well. Taking advan- 
tage of the sixfold rotational symmetry of the lattice, we 
need to consider only a small subset of spins (marked cri 
through cTg by filled dots in Fig. 1) to provide the mag- 
netization of the whole particle. Each of the spin sites, 
which stand for the atomic moments in the nanoparticle, 
are described by Ising spin variables that take on values 
ai = ±1. The nanoparticle is assumed to be very small, 
so that it constitutes a monodomain. In addition, accord- 
ing to experimental results ^J, an uniaxial anisotropy is 
implicit in the model, with the easy axis coinciding with 
the Ising spin direction. Therefore the energy of the sys- 
tem is written as 



H = — ^ Jij(7i(Tj — H^y^aj 



(1) 



<«j> 



where only the exchange coupling between nearest neigh- 
bors is considered, with as the exchange coupling con- 
stant, and H is the applied magnetic field. 

Having defined the relevant spins that describe the sys- 
tem we start by writing the dynamical equations for the 
average spin magnetizations 



a 

and the pair correlation functions 

Tijit) OiOj ^ ajajPia, t), 



(2) 



(3) 



where the index i, as well as j, designates a spin site and 
the sum over a accounts for all the spin configurations. 
Since we restrict our model to take into account only pair 
correlations, the probability of a given spin configuration 
is written as 



[ i i<J ) 

where iV is the total number of spins. 

The time evolution of the system is obtained by means 
of the master equation for the probabilities 



with w{a',a) as the transition probability per unit time 
from the spin configuration a' to configuration a. Here 
we assume that the dynamics of the system is governed 
by single-spin-fiip transitions, which can be mathemati- 
cally described by the Glauber formalism fiq] 



W^{cr^) = — {1 - (7,7} 



(6) 



with 



y( J, /3, H) = tanh [(^ J.^a, + H)/kBT] (7) 



= tanh (Ai), 



(8) 



where fc^ is the Boltzmann constant, T is the absolute 
temperature and r is the relaxation time for a single spin. 
This dynamical rule is chosen so that the transition prob- 
abilities WiijJi) depend on the spin value of Ui as well as 
on the values of its nearest neighbors. 

Using Eqs. (|^) to (^ we obtain after some alge- 
braic manipulation the dynamical equations for mi(t) 
and Vij (t) 



dmi{t) mi(t) 1 

dt T T 



^tanh(A,)P,(a,t) , 



dt 



(9) 



(10) 



- ^ [ai tanh (A^-) + aj tanh (A,;)] {a, t) 



In principle, our task is to integrate these equations to 
obtain the time evolution of the variables of interest. 
Within the Bethe-Peierls approximation we assume 
that only the correlations between pairs of nearest neigh- 
bors are nonzero and write the following equations for 



Ua'Pi'^^,<^o 



(11) 



where Ni is the coordination number of the lattice site i 
and for (cr) 



P.A<^) 



X ]^P(o-.i,cr„) J|P(o-,,cr/3) 



(12) 
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It is important to notice that Ni is different for the sur- 
face spins. P{<Ji) and P(ai,crj) are obtained from ex- 
pression (Q) taking as nonzero only the nearest neighbor 
correlation functions 



(13) 



P{ai,aj) = -{I + rtiitJi + m.jaj + Tijaiffj) . (14) 

Substituting expressions (^Tj) and into the dy- 
namical Eqs. (||) and (|lo|) yields, after a lengthy but 
straightforward manipulation, a closed system of equa- 
tions for the spin magnetizations {mi{t)) and nearest 
neighbor pair correlation functions (rij{t)) that can be 
solved numerically 

dmAt) , , 



dt 



k=0 



^ ^ '^V'^Z^m}) + %u'^W,i{N,}) } , (15) 

'^i Vi 



(16) 



with 



Af,-1 



V^Z.iiN,}) 



k=0 



fc=0 



(17) 



and obtained by the simple exchange of the indices 
i and j in Equation (|l^). Refer to the appendix for a 
definition of the symbols. We also define the products 



N, 



Z,{{N,}) = l[[z^], 

a. 



upon which act the operators 
V 



^Zai]j ■ ■ ■ [Zak]j 



Ql<...<Qfc 



E 



(18) 
(19) 

(20) 
(21) 
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that generate all the possible configurations of the clus- 
ter. In the case of Eq. (^ the independent indices ak 
stand for all the nearest neighbors (1 < afc < 6) of a 
given spin aj, as depicted in Fig. 2(a). When calculat- 
ing the pair correlation = {ai<jj), to avoid counting 
the pair twice 1 < afe < 5, as shown in Fig. 2(b) for 



spins o-i and Uj. For A; = 0, V° = 1 and = 1. The 
parameter is related to the transition rate of a spin 



tk = tanh 



f2J{k-3) + H 



ksT 



(22) 



The results comprehending Eqs. (y_5|) throughout ( |21| ) 
correspond to the case Jy = J. To account for the more 
complex case of different in heterogeneous nanopar- 
ticles V*^ and U'^ have to be modified to incorporate a 
generalized parameter tk- This is done in the appendix. 



III. RESULTS AND DISCUSSIONS 



In this section we describe the properties of the model 
and show that it accounts for many of the phenomena 
presented by magnetic nanoparticle systems. 

We start by discussing the properties of the hysteresis 
of homogeneous ferromagnetic nanoparticles, as yielded 
by the present model. For this purpose consider Fig- 
ure 3(a) where it is shown hysteresis loops for 3 different 
nanoparticle sizes i? = 1 , 2 and 11. i? is the num- 
ber of shells and can be associated with the radius of 
the nanoparticle. We define the total magnetization of 
the nanoparticle M = mi(t)/N, where the sum runs 
over all the spins of the nanoparticle (iV). The small- 
est nanoparticle (i? = 1) has vanishing coercive field and 
its magnetization is well fitted by a Langevin function, 
whereas the bigger particles (i? = 2 through 11) exhibit 
typical hysteresis loops of ferromagnetic nanoparticles. 
Figure 3(b) shows the dependency of the coercive field 
{He) on the number of shells of the nanoparticle, by plot- 
ting He as a function of R^^-^, for the reduced tempera- 
ture ksT = 3J. This graph demonstrates that the fimc- 
tional relation He ^ R^^'^ fits very well the results of our 
calculations for nanoparticles of size 2 < i? < 10, with 
the exception of the R = 1 and i? = 11 cases. The for- 
mer is in the superparamagnetic regime [l^ , as indicated 
by Fig. 3(a), whereas the last already starts to behave 
as an infinite system. The same behavior exhibited by 
He in Fig 3(b) is obtained for the reduced temperature 
ksT = 1.5J. According to our calculations, for i? > 10 
the Curie temperature (Tc) of the nanoparticle coincides 
with the Te of an infinite cubic system, as yielded by the 
Bethe-Peierls relation coth [J/fc^Tc] = 5 [||. Calcula- 
tions for a square lattice showed that He ^ L^^-^ for 
small clusters, where L stands for its lateral dimension. 
Finally, we point out that the mechanism of magnetiza- 
tion reversal that takes place in our nanoparticle is not 
entirely coherent. Because of their smaller coordination 
number, the average magnetization (m^) of the outside 
spin shells decrease by the action of the reverse field H, 
consequently also decreasing the magnetization of their 
neighboring shells and creating a nonuniform magnetiza- 
tion profile. Eventually, when H reaches He all the spins 
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flip together. This process should be responsible for eas- 
ier magnetization rotation and lower switching fields than 
those given by the Stoner-Wohlfarth model. 

In the remainder of this section we consider the proper- 
ties of heterogeneous nanoparticles, which are more com- 
plex systems consisting of ferromagnetic and antiferro- 
magnetic (or ferrimagnetic) regions coupled with each 
other. Examples of such systems are oxide-coated cobalt 
(Co) or iron (Fe) nanoparticles among many others 
that exhibit the exchange anisotropy effect . When 

cooled in the presence of an external magnetic field the 
soft ferromagnetic core of the particle aligns with the ap- 
plied field but its outer antiferromagnetic (AF) surface 
gets ordered only when the temperature of the system is 
lower than Tn (Neel temperature) The coupling of 
the ferromagnetic core with the antiferromagnetic surface 
therefore produces a unidirectional exchange anisotropy 
that shifts the hysteresis loop to higher or lower mag- 
netic fields, depending on their mutual orientation. To 
obtain this effect we consider an 8-shell nanoparticle with 
the following characteristics: the exchange interactions in 
the core (i? = to 7) are given by the coupling parameter 
Jc = 1, the outmost shell {R = 8) is also ferromagnetic 
with an exchange constant Js = 5 and the coupling be- 
tween the R = 7 and R = 8 shells is due to an antiferro- 
magnetic exchange constant Jc-S = —4. The hysteresis 
we obtain with our formalism is shown in Figure 4(a), 
which presents a clear exchange shift. Figure 4(b) shows 
the magnetization curves of individual spins within the 
nanoparticle: in the center of the nanoparticle (Sq), in 
the sixth shell (5*6), in the seventh shell (Sr) and in the 
surface shell (S's). As evidenced by this graph, the sur- 
face shell and the layer just below it are locked in a very 
stable AF configuration, whereas the average magneti- 
zation of the inner ferromagnetic shells vary gradually 
towards the center of the particle. For the sake of com- 
parison, we have also calculated the hysteresis for the 
same system, but using a simple mean field approach, 
whose results are shown in Figures 4(c) and 4(d). For 
the standard mean field equations each spin experiences 
only an effective field that can be seen as an average in- 
teraction with all the other spins and, consequently, the 
exchange shift disappears. Looking at the individual spin 
magnetizations (Fig. 4(d)) it is evident that there is no 
correlation between the surface and core domains. 

Another important experimental technique used to in- 
vestigate the properties of magnetic nanoparticles is the 
field cooling (FC) and zero field cooling (ZFC) magneti- 
zation measurements. For an ensemble of non-interacting 
homogeneous particles the blocking temperature (Tb) 
is related to the magnetic anisotropy constant and the 
volume of the nanoparticle J 1^ . For a system of in- 
teracting nanoparticles see p]. The blocking tempera- 
ture, where the maximum of the ZFC curve occurs, in- 
dicates the point in which the thermal energy is compa- 
rable to the average anisotropy barrier in the nanopar- 
ticles, and beyond it the magnetization decreases with 
temperature. In the following we show a deblocking 



behavior that is caused only by the exchange interac- 
tion between antiferromagnetic and ferromagnetic re- 
gions in a single nanoparticle. For that purpose consider 
a 6-shell nanoparticle whose core is predominantly AF, 
Jc = —0.5 Jo for i? = to 5, with a ferromagnetic surface, 
J5 = 2Jo for i? = 6, which is ferromagnetically coupled 
with the core, Jcs = ^Jo between R = 5 and i?=6. The 
energy parameter is Jq = 1 and the time taken to sweep 
over the whole temperature interval is ts = lO^r, with 
T appearing in Eqs (||) and (|o|). The calculations were 
performed in accordance with the standard FC-ZFC ex- 
perimental practice The obtained FC (solid line) 
and ZFC (dashed line) magnetization curves are shown 
in figures 5(a) through 5(f) as a function of the reduced 
temperature, for values of the applied field that vary from 
H = 0.05|Jo| to H ^ 0.0001|Jo|. M is the normalized 
magnetization of the entire nanoparticle, as previously 
defined. Because it corresponds to a single particle, Tb 
is sharply marked and does not exhibit the characteristic 
dispersion caused by the size distribution of the parti- 
cles. According to our calculations, Tb is not strongly 
dependent on Jc however it does depend on Jsc and Jc- 
Figure 5 shows that the blocking temperature occurs just 
about the range Jcs ^ ksTB < Js, which is confirmed 
for different nanoparticle parameters (size, exchange con- 
stant and ts). Despite being observed to take place for 
nanoparticles of different sizes, the blocking temperature 
effect is very sensitive to the relative strengths between 
the ferromagnetic and antiferromagnetic exchange con- 
stants. For clarity. Figure 6 displays the relation between 
the average magnitude of the energies associated with the 
parallel spin alignment {Up + Ur) and anti-parallel spin 
alignment {Uaf) configurations within the particle. 



C/f = ( ^ Js O-iCTj + ^ Jcs CfkCJl) 

<i,j> <kd> 

= Js ^ Tij + Jcs ^ rti , 



<ij> 



<k,l> 



■i i 

Uaf = Jc ct^ctj) = Jc ^ 



(23) 
(24) 
(25) 



<ij> 



<«j> 



in terms of which we define the differential exchange en- 
ergy Axe = {Up + Uh- Uaf)/{Uf + Uh + Uaf)- Solid 
and dashed lines represent FC and ZFC processes, re- 
spectively, and the arrows indicate their dynamics. The 
contribution of the Uh term to the total energy of the 
nanoparticle is very small, however it is responsible for 
the deblocking behavior at Tb- At the high tempera- 
ture regime the average exchange energy is nonzero, be- 
cause of a remaining finite short range correlation be- 
tween nearest neighbor pairs. Most of the energy is due 
to the ferromagnetic ordering. During the ZFC (dashed 
line) the AF core overwhelms the ferromagnetic surface 
and forces the magnetization of the whole nanoparticle 
to vanish. Then, in the subsequent heating, even a very 
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small applied field is enough to trigger the alignment of 
the ferromagnetic surface, which is evidenced by the rel- 
ative increase of the ferromagnetic energy. On the other 
hand, during the FC process the applied field is able to 
magnetize the nanoparticle at higher temperatures, be- 
fore the AF correlations of the core become too strong. 
Thus, the interplay between AF and F exchange ener- 
gies gives rise to the hysteresis in temperature that we 
observe in the FC-ZFC calculations. 



(30) 



where [zi]j = [zj]i and [wi]j = [wj]i. As already men- 
tioned in the main text, the parameter tk in Eqs. ([l5| ) 
and (nq) is related to the transition rate of a spin 



tk = tanh — — 



(31) 



IV. CONCLUSIONS 

We have developed a new approach to calculate the 
classical magnetic properties of single-domain nanoparti- 
cles. Treating the nearest neighbor pair correlation func- 
tions explicitly we were able to obtain a simple set of 
equations for the magnetization at any point in the par- 
ticle. In addition, Eqs. (|l^) and ( p^ ) can also account for 
the inclusion of interactions between next nearest neigh- 
bors in a straightforward manner, by incorporating the 
corresponding Jij into the operators (|3^ ) and (|33| ) and 
increasing the coordination number Nj of each spin. The 
method can be applied to different lattice geometries and 
extended to 3D systems. The formalism treats correctly 
the short-range correlations that arise in heterogeneous 
nanoparticles, which have regions of different magnetic 
properties, e.g., oxidized ferromagnetic particles. Such 
spin correlations were evidenced by the observation of 
the exchange shift in hysteresis loops and by the blocking- 
deblocking effect in the FC-ZFC calculations. 
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for Jij = J and in accordance with the definitions given 
in the main text. For the case of general exchange con- 
stants (Jij), the term tk must be generalized and incor- 
porated into the operators V'' and 



2-^ 1 



tanh 



E 



Q;i<...<Q!i; 



tanh 



so that, for k — Q 



V" = tanh 



W = tanh 



knT 




(32) 



(33) 



(34) 
(35) 



VI. APPENDIX 



Here we define some of the notation used in section ||. 
Based on Eq. (|), equations (|l|) and (|l|) are written 
out in terms of the following simple functions 



1 - rrii 



1 + TOi + Urij + Ti 
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I - mi 


- ruj -t- Ti-j 
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I - nil 


+ rrij - Tij 



(26) 



(27) 



(28) 



(29) 



[1] E.G. Stoner and E.P. Wohlfarth, Philos. Trans. R. Soc. 

London, Ser. A 240, 599 (1948); reprinted in IEEE Trans. 

Magn. 27, 3475 (1991). 
[2] S.A. Makhlouf, F.T. Parker, F.E. Spada, and A.E. 

Berkowitz, J. Appl. Phys. 81, 5561 (1997); R.H. Ko- 

dama, S.A. Makhlouf, and A.E. Berkowitz, Phys. Rev. 

Lett. 79, 1393 (1997). 
[3] W.H. Meiklejohn and CP Bean, Phys. Rev. 105, 904 

(1957); S. Banerjee, S. Roy, J.W. Chen, and C. Chakra- 

vorty, J. Mag. Mat. 219, 45 (2000). 
[4] R.H. Kodama, A.E. Berkowitz, E.J. McNifT Jr, and S. 

Foner, Phys. Rev. Lett. 77, 394 (1996); B. Martinez, X. 

Obradors, LI. Balcells, A. Rouanet, and C. Monty, Phys. 

Rev. Lett. 80, 181 (1997). 



5 



[5] R. Prozorov, Y. Yeshurun, T. Prozorov, and A. 

Gcdankon, Phys. Rev. B 59, 6956 (1999). 
[6] W. Wernsdorfcr, E.B. Orozco, K. Hassolbach, A. Bonoit, 

D. Mailly, O. Kubo, H. Nakano, and B. Barbara, Phys. 

Rev. Lett. 79, 4014 (1997); E. Bonet, W. Wernsdorfer, 

B. Barbara, A. Bcnoit, D. Mailly, and A. Thiaville, Phys. 

Rev. Lett. 83, 4188 (1999). 
[7] M. Jamet, W. Wernsdorfer, C. Thirion, D. Mailly, V. 

Dupuis, P. Melinon, and A. Perez, Phys. Rev. Lett. 86, 

4676 (2001). 

[8] J. Shen, R. Skomski, M. Klaua, H. Jenniches, S.S. 
Manoharan, and J. Kirschner, Phys. Rev. B 56, 2340 
(1997). 



[9] M. Plischke and B. Borgorson, Equilibrium Statistical 
Physics, 2nd edition (World Scientific, Singapore, 1994). 

[10] R.J. Glauber, J. Math. Phys. 4, 294 (1963). 

[11] K. Huang, Statistical Physics, 2nd edition (Wiley, New 
York, 1987). 

[12] B.D. Culht, Introduction to Magnetic Materials 

(Addison- Wesley,London, 1972). 
[13] W.H. Meiklejohn, J. Appl. Phys. 33, 1328 (1962). 
[14] R.C. O'Handley, Modern Magnetic Materials: Principles 

and Applications (Wiley & Sons, New York, 1999). 
[15] R. Sappey, E. Vincent, N. Hadacek, F. Chaput, J. P. 

Boilot, and D. Zins, Phys. Rev. B 56, 14551 (1997). 



FIG. 1. Schematic representation of an hexagonal nanoparticle exhibiting 3 shells of spins (small circles). The representative 
spins that compose the set of equations for this particle are represented by black circles. The dashed lines evidence the sixfold 
symmetry. 

FIG. 2. (a) A cluster formed by spin aj and its 6 nearest neighbors, (b) A cluster formed by the pair aiaj and their nearest 
neighbors. 

FIG. 3. (a) Hysteresis loops for various nanoparticle sizes: 7i = 1,2 and 11. The reduced temperature is T = SJ/ks, where 
J = 1 is the exchange coupling constant and kn is the Boltzmann constant, (b) The coercive field {He in units of J) plotted 
as a function of R~^'^. From right to left, the dots correspond to ii = 2, 3, 11 and the straight line results from a lineax fit 
to the data. 

FIG. 4. (a) Hysteresis loop for a nanoparticle of radius R = 8. Jc; = 1 for R < 7, Js = 5 in the surface R — 8, and 
Jc-s = —4 couples the core and surface regions, (b) Individual magnetization curves for the center spin (^o), a spin in the 6th 
shell (Se), a spin in the 7th shell (Sr) and a spin in the surface shell (Sg). (c) Same parameters as in Figure (a) using a simple 
mean field approach, the exchange shift disappears, (d) Same as Figure (b) for simple mean field. The reduced temperature is 
ksT = 3.5 Jc for all curves. 

FIG. 5. FC (solid) and ZFC (dashed) curves of the total normalized magnetization (M) as a function of the reduced temper- 
ature fcsT/Jo. Different values for the applied magnetic field are considered, from H = 0.05|J()| in (a) through H — O.OOOljJoj 
in (f), as indicated in the figures. The nanoparticle consists of R = 6 shells, with an AF core (Jc = —0.5 Jo, for < J? < 5), 
a ferromagnetic surface (Js = 2Jo within R = 6) and a ferromagnetic coupling between core and surface {Jcs = Uo, between 
R = 5andR = 6). 

FIG. 6. Differential exchange energy Axe = {Uf + Uh — Uaf)/{Uf + Uh + Uaf) as a function of the reduced temperature 
ksT/Jo. The applied field is H = 0.001] Jo]. 
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FIGURE 1 
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FIGURE 2 
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FIGURE 4 
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FIGURE 6 
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